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Abstract 

The microscopic effective pairing interaction in the 1 S'o-channel is investigated 
for two different values of the chemical potential jjl starting from the separable form 
of the Paris NN-potential. It is shown that, within a high accuracy, this effective 
interaction can be approximated by the off-shell free T-matrix taken at the negative 
energy E = 2\i. 



Recently |1[] the microscopic effective pairing interaction Vf ff in the 1 S'o -channel 
was found by solving the Bethe-Goldstone equation for semi-infinite nuclear matter 
without any form of local approximation. The separable representation [0, || of the 
Paris potential Q was used which simplifies the problem significantly and makes 
it possible such direct solution for non-uniform systems with the use of the mixed 
coordinate-momentum representation. However, the procedure turned out to be very 
cumbersome and the effective interaction obtained has rather complicated form not 
convenient for applications. In this paper we show that V^ s , with a good accuracy, 
can be approximated by the free T-matrix taken at the energy E = 2//. This 
quantity is much simpler than V^ ff and can be easily found for a fixed value of \i. 

The many-body theory form of the gap equation |5|, |G|] 

A = VGG S A (1) 

is used in which explicitly takes into consideration the particle-particle propagator 
A s = GG S in the superfluid system. Within the Bethe-Brueckner approach, the 
irreducible particle-particle interaction block V is approximated by the free NN- 
potential. The effective interaction is associated with the splitting of the complete 
Hilbert space S into two domains, the model subspace S and the complementary 
one S'. As a result, the two-particle propagator is represented as a sum A s = A S Q +A'. 
It is supposed that the superfluid effects can be neglected in the S"-subspace and 
therefore the upperscript "s" is omitted in the second term. 

The gap equation (1) can be rewritten in the model subspace: 

A = V^A, (2) 

where V^ ff obeys the following equation: 

V c p ff = V + VA'V» S . (3) 

In H the model space So is defined in such a way that it involves all the two- 
particle states (A, A') with the single-particle energies e\,ey both negative. In this 
case, the complementary subspace 5" involves all the two-particle states with positive 
energies E\,E\i and two-particle states with one energy positive and the second one 
negative (but greater than //). 

The separable form |2|, £| of the Paris potential is as follows: 

V(k,k>)=J2KMk 2 Mk' 2 ), (4) 

ij 

where form factors gi{k 2 ), (i = 1,2,3) are rational functions. This choice was 
adopted earlier in the Brueckner-type calculations for infinite nuclear matter 0, || . 



It should be noted that the original normalization [fj, [| of the expansion (4) was 
changed in [T] in such a way that the identity (?j(0) = 1 holds true. Then the 
absolute values of the Ay-coefficients give direct information on the strength of 
the corresponding terms of the force. Their values (in MeV-fm 3 ) are as follows: 
A n = - 3659, Ai2=2169, A 22 = - 1485 and Ai 3 = - 23.6, A 23 =57.6, A 33 =17.2. As it 
is seen, the strengths of all the components containing only the indices i = 1, 2 are 
much stronger than those with the index i = 3. Therefore the latter will not be 
considered for a qualitative analysis. Of course, in the calculations all the terms A^ 
are considered. 

The separable representation of the iViV-potential leads to a similar form of V^ ff 
which, in the notation of |l| is as follows: 

V e p ff (£rj_, k'l; x 1 ,x 2 , x 3 , x 4 ; E) = J2 Kj(X 12 , X u ; E)gi(k\, x 12 )gj(k'l, x u ). (5) 

ij 

Here the center-of-mass and relative coordinates in the x-direction are introduced 
(Xl2 = (x\+x 2 )/2, x\ 2 = x\—x 2 , etc.), and g%(k±,x) stands for the inverse Fourier 
transform of the form factors gi(k± + k%) in the x-direction. 
The coefficients A^ obey the set of integral equations: 

Ay(Xi2, X 3i ] E) = Xij6(X 12 — X 34 )+ 

+ 2^ A,/ / dX 5e Bi m (Xi 2 , X 56 ; E) A m j(X 56 , X 34 ; E), (6) 

lm 

where B\ m are given by integrals of the propagator A' with two form factors. Their 
explicit form is as follows: 



u (v v t?\ S^' f d 2 ^± G nn , (kj_,X 12 ) G'n* n (kj_,Xu) 

^ J ( 27r ) 2 E-S X - Sy 

G l nn ,(k 2 ± ,X 12 ) = / dx 12 gi(k 2 ± ,x 12 ) y n (X 12 +x 12 /2)y n/ (X 12 -x 12 /2), 



(7) 



where A=(n, kj_), e\=e n + k 2 _/2m, e n , y n stand for the energies and wave functions 
of the one-dimensional Schrodinger equation, respectively. The prime in the sum of 
eq. (7) means that the summation is carried out over (A, A') which are not included 
in the model space. It is convenient to put in this sum e\<ey (multiplying the result 
by the factor 2). Then the sum contains all the states with e\>fi, ey>0. 

To overcome a problem of a slow convergence of the integrals of eq. (7), a renor- 
malization can be made in terms of the free T-matrix: 

T = V + VA°T, (9) 

where A is the propagator of two free particles. 



Eqs. (6), (7) hold true for the T-matrix (with the substitution T — > A, B° — > 
B), but now the coefficients Tij(X, X'; E), B^(X, X'\ E) of the separable expansion 
depend only on the difference t — X — X ' of the CM coordinates: 

T tJ (t; E) = \ 5 8{t) + £ A a / dt ' B imit ~ t'; E) A mj (t'; E), (10) 

Ira 

The renormalized equation for the effective interaction in a compact form reads: 
Aij = Tij + Y^ T ii (Bim - B lm ) T mj (11) 

lin 

The kernel of this equation converges much faster than the original one. Of 
course, the problem of the slow convergence at large momenta does not disappear. 
It passes to eq. (10), but in this case the difficulty can be overcome much more 
easily. Indeed, we are dealing now with the one-dimensional vector Tij(t) instead of 
the two-dimensional matrix Aij(X,X'). It is convenient first to find the T-matrix 
in the momentum representation by solving the following set of equations: 

T i3 {P x - E) = \ %j + £ \ u B? m (P x ; E) T mj (P x ; E), (12) 

lm 

where 

B lm {P x , E)-J ( ^ )3 - _ p2jAm _ k2/m [16) 

Then Ty(£) can be found from Tij(P x ) with the help of the inverse Fourier transfor- 
mation: 

7 dp 

Tij{t; E)= J -^ T l3 {P x - E) eM-iPxt). (14) 

— oo 

The form factors Qi in eq. (13) are rational functions on k 2 0,0]) namely, 
combinations of the Yukawa function and their derivatives with different masses 
Pin (n — 1, - - • , 4). This integral can be evaluated analytically, but such calculation 
is very cumbersome because of a huge number (~70) of particular terms appearing 
in the integrand. We prefer, following [|l|], to integrate it numerically, with the cut-off 
momentum k c = 60 mi" 1 which guarantees an accuracy better than 1%. 

The Fourier integral (14), after separating the constant term Ay, was calculated 
in |1| by direct integration along the real P x -axis. This method works well at small t, 
but for t > (2 -=- 3) fm the integrand contains rapidly oscillating factors multiplied by 
a slowly falling function (Tij(P 2 ) — Ay) — l/P x , which makes the convergence very 
poor. To obtain a reasonable accuracy at t=(4 -h 5) fm, the large value of the cut-off 
momentum P%. = 3000 fm^ 1 was taken in |l| with very small integration step. It is 
much better, in accordance with the general recipe of integrating rapidly oscillating 



functions, to integrate eq. (14) in the complex plane of P x . The integration contour 
C can be closed in the upper half-plane and deformed to a form convenient for 
numerical integrating (Fig. 1). As it can be readily shown, in the case of rational 
form factors gi all the singularities of Tij(P x ; E) are the poles (simple and multiple) 
located on the imaginary axis (symmetrically to the origin). Their position P^—ij a 
is defined by various combinations of the masses (3 in and the /x-dependent parameter 
7o = •>/— 8m/x. The contour C surrounds all these poles, each one yielding a falling 
exponent exp(— 7 Q |t|). Besides, the procedure of numerical integration of eq. (13) 
produces some additional "false" poles depending on k c . Though their contribution 
is negligible, the contour surrounds them also for sake of consistency. In practice, 
we used the contour C with values of parameters a = 2fm _1 , b = 120 fm -1 . In this 
case it is placed far enough from all the poles and the integral can be calculated 
without numerical problems yielding the correct result for any value of t. 

We describe the procedure of calculating the off-shell T-matrix in detail because, 
as will be shown, the difference B — Bq in eq. (11) turns out to be to be rather small 
and the solution of this equation with a good accuracy coincides with the T-matrix. 
Therefore the evaluation of this quantity is of primary importance. 

In this paper, we go from the semi-infinite system to a more realistic geometry 
of a finite slab within the Saxon- Woods potential well. The parameters are chosen 
in such a way to reproduce qualitatively those of heavy nuclei in the lead region: 
the width of the slab is 2L = 16 fm, the depth of the well Vq = —50 MeV and the 
diffuseness parameter d = 0.65 fm. In the slab case, all the above equations are valid 
with substitution of the slab wave functions and energies in eqs. (7), (8). Instead of 
direct solution of these equations, the Local Potential Approximation (LPA) is used 
which was previously proved to be a reliable approximation for semi-infinite nuclear 
matter [3J. Within the LPA, the exact values of Bi m (Xi, Xi; E) are replaced by 
the set of those Bf^(t, E; V[X]) for infinite nuclear matter put in the potential well 
V[X], Here X = (X\ + X 2 )/2 is the average value of the two CM coordinates. As 
it is shown in [Q, comparing exact solution for semi- infinite matter with the LPA 
prediction, the latter works with an accuracy of few percent even in the surface 
region. It is natural to suppose that in the slab case the LPA works also sufficiently 
well. The LPA procedure is as follows. At a fixed value of the chemical potential /^, 
it is necessary, first, to calculate the set of vectors B^(t, E = 2/i; Vi) (Vi=5V-(i — 1)), 
and then, for every value of (X,t), to find Bj^ A (X 1 ,X 2 ) by interpolating values of 
B'l^(t; V) with nearest to V(X) values of V. 

Let us first take /i= - 8 MeV. Shown in Fig. 2 is the difference (B}£(V) - Bf m ) 
for three main components Im, in comparison with the free propagators Bf m , for 
the case of the maximum value of the potential depth V = 50 MeV. Of course, 
for smaller values of V this difference is even less. To analyse the origin of this 
smallness, let us separate the sum of eq. (7) into two parts, the first one with both 



positive energies and the second one with one negative energy, and compare each 
of them with the corresponding sum for the free propagator. It should be stressed 
that within the LPA eq. (7) is very similar to that for Bi m because the plane waves 
in the constant potential Vq=V(X) stand for y n in this case. Every term of the first 
sum corresponds directly to the one of Bf m with the same energy denominator. It 
can be easily seen that the numerator of the in-matter term is less than that of the 
free one because the matrix elements of the form factors are calculated for bigger 
values of momenta (q = \/p 2 + 2m V in the in-matter case in comparison with p 
in the free one). Since all the form factors gi(p 2 ) fall with p, this part of B}^(V) 
(in absolute value) is less than Bf m (at 20 -j- 30 %). The second part of the sum 
under consideration has no analogue in the free propagator. Though it contains an 
essentially smaller phase space that the first one, its value turns out to be significant 
due to small value of the energy denominators. Its contribution to the absolute value 
of Bjm(V) is even larger than the difference of the first term of Bi m and Bf m , so the 
absolute value of Bi m exceeds that of Bf m , but only a bit. The size of the difference 
is about 10%. 

Fig. 3 shows similar comparison of the effective interaction V^ [V=50 MeV] and 
the free T-matrix. Of course, the <5-terms of both amplitudes (see eqs. (6), (10)) are 
extracted. For more descriptive comparison of the "exact" LPA effective interaction 
and the free T-matrix we display in Fig. 4 their zero-order moments: 

A^X) = JdtA tJ (X-t/2,X+t/2), (15) 

and the same for T^-. As it is seen, the difference between these averages is again very 
small for the 11-component, but for other components it is greater, namely ~ 20% 
for the 12-component and ~ 30% for the 22-one. However, these two deviations are 
of the opposite sign and, as we show just below, almost compensate each other in 
the sum at values of k 2 , k' 2 which are important for the pairing problem. 

It can be demonstrated by analysing the localized form of the effective interaction 
which was introduced in [|TJ: 

V&P0 = Eh(X)9i(kUX)) gj (kUX)), (16) 

ij 

where kp(X)=2m([/,—U(X)) if //— U(X) > and kp(X)=0 in the opposite case. In 
f|], it is argued that, for states in the Fermi surface vicinity, the effective interaction 
can be approximately replaced by V^(X). It is drawn in Fig. 5, together with T F 
which is defined in a similar way. One sees that the compensation discussed above, 
indeed, occurs and the difference between the two curves is practically negligible. 

To imitate the situation in the vicinity of drip-lines, all the calculations were 
repeated for /i= — 4MeV. Results are completely similar to those for jjl— — 8MeV. 



They are shown in Fig. 6 for the B\\ component and in Fig. 5, for the local form 
of the effective interaction V^ S (X). Thus, this simple approximation of Vg ff with 
the free T-matrix can be used also for predicting properties of nuclei nearby the 
drip-lines. It should be mentioned that one approximation used in this analysis 
looks doubtful for small values of fi. This is neglecting all the pairing effects in the 
complementary space. This drawback can be improved by small change of the model 
space by addition to it a set of states with small positive energies E\ < E where 
Eq is of the order of several MeV. Estimates show that in this case the difference 
between V^ s and the T-matrix becomes even smaller. 

Thus, the free off-shell T-matrix taken at the energy E = 2/i, indeed, is a good 
approximation for the microscopic effective pairing interaction. For the separable 
representation of the Paris iViV-potential it can be easily found by solving eqs. (11)- 
(13). As far as we deal with comparatively small shift from the mass shell, all the 
realistic iViV-potentials must give approximately the same predictions for the T- 
matrix. 
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Figure 1: The integration contour for the inverse Fourier transformation (14). 
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Figure 2: The free propagators —Bf m (solid lines) and the differences —(B^[V = 
50 MeV)] —Bf m ) (dotted lines) calculated for \i = — 8 MeV versus relative coordinate 
t = X-X'. 
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Figure 3: The /m-components of the effective interaction V^ [V = 50 MeV] (solid 
lines) and the free T-matrix (dotted lines) calculated for \x = —8 MeV versus t. Sign 
is changed for i,j — 1, 1 and 2, 2. 
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Figure 4: The zero moments Aij(X) (solid lines) and the values Tij (dotted lines) 
for fi = -8 MeV. 
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Figure 5: The effective interaction Vj^(X) (solid lines) and T?(X) (dotted lines) 
for two values of the chemical potential \x (the number corresponds to \\x\ in MeV). 



13 



3 i\/r^/-i fv^-4 



B n (t) ■ 10 3 , MeV" 1 - fm 




1 ' 



2.0 t, fm 



Figure 6: The same as in Fig. 2 for I, m — 1, 1 calculated for // = —4 MeV. 
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